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Abstract 

Nonlinear variants of quantum mechanics can solve tasks that are impossible in standard 
quantum theory, such as perfectly distinguishing nonorthogonal states. Here we derive the 
optimal protocol for distinguishing two states of a qubit using the Gross-Pitaevskii equation, a 
model of nonlinear quantum mechanics that arises as an effective description of Bose-Einstein 
condensates. Using this protocol, we present an algorithm for unstrnctured search in the Gross- 
Pitaevskii model, obtaining an exponential improvement over a previous algorithm of Meyer 
and Wong. This result establishes a limitation on the effectiveness of the Gross-Pitaevskii 
approximation. More generally, we demonstrate similar behavior under a family of related 
nonlinearities, giving evidence that the ability to quickly discriminate nonorthogonal states and 
thereby solve unstructured search is a generic feature of nonlinear quantum mechanics. 

1 Introduction 

Linearity is an essential feature of quantum mechanics whose violation can have dramatic oper¬ 
ational consequences. In particular, Abrams and Lloyd [3] showed that in a model of nonlinear 
quantum mechanics due to Weinberg [11], one can exponentially increase the angle between quan¬ 
tum states, distinguishing states separated by an angle e in time 0(log J). 

This ability has strong computational implications, including a fast algorithm for unstructured 
search (and thus a fast algorithm for any problem in NP). In the unstructured search problem, one 
aims to find a marked item using a black box that determines whether any given item (out of N 
possible items) is marked. Using the ability to distinguish nonorthogonal states in the Weinberg 
theory, Abrams and Lloyd gave an algorithm for unstructured search that uses time 0(log A) and 
only a single query. 

More recent work of Meyer and Wong [9] considered the computational power of the Gross- 
Pitaevskii equation, a nonlinear Schrodinger equation that provides an effective description of 
Bose-Einstein condensates. Specifically, they considered a discrete version of the Gross-Pitaevskii 
equation. 
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for a unit vector |?/)) G span{|x)}, where the coefficient g quantifies the strength of the nonlinearity. 
They presented an algorithm for unstructured search that simultaneously applies an oracle Hamil¬ 
tonian, an input-independent driving Hamiltonian, and the Gross-Pitaevskii nonlinearity. Their 
algorithm uses time 0(min{so it can (unsurprisingly) solve the search problem arbi¬ 
trarily fast using an arbitrarily strong nonlinearity. However, by considering the resources needed 
to measure the system at a precisely chosen time, they argued that a reasonable implementation 
would have complexity giving only a modest improvement over the 0{-\fN) complexity 

of Grover’s algorithm [7]. 

In this paper, we apply the original approach of Abrams and Lloyd to the Gross-Pitaevskii 
nonlinearity. In Section 2, we find the optimal protocol for distinguishing two states of a qubit by 
adding a carefully-chosen driving term. The optimal procedure distinguishes states with overlap 

1 — e in time O(Mog^), recovering essentially the same behavior as the Weinberg model. In 
Section 3, we apply this result to the unstructured search problem. We present an algorithm with 
complexity 0(min{^ \og[gN), \/iV}), exponentially improving the upper bound of Meyer and Wong. 

Furthermore, we give a lower bound of H(min{^, showing that our algorithm is optimal up 

to a logarithmic factor. 

It would be interesting to find a variant of quantum theory that allows a polynomial speedup over 
quantum computation, but no more (see [2] for some recent progress in this direction). However, 
our results suggest that the Gross-Pitaevskii model does not provide such an example. We also 
consider a broad class of nonlinearities that generalize the Gross-Pitaevskii model and show that 
all such models lead to essentially the same behavior. This supports the view that the ability to 
exponentially increase the angle between states is a generic feature of nonlinear quantum mechanics, 
as previously suggested by Abrams and Lloyd [3] and Aaronson [I] . 

In light of the dramatic consequences for information processing, it seems unlikely that quantum 
mechanics is fundamentally nonlinear. However, information-processing protocols using the Gross- 
Pitaevskii model can be used to place limitations on the validity of the underlying approximation 
[9]. We discuss this point further in Section 4, where we use our state discrimination protocol to 
show that for an AA-atom Bose-Einstein condensate, the Gross-Pitaevskii approximation can only 
hold up to time 0{^ logAA). 

Finally, we conclude with an open problem regarding higher-dimensional nonlinear state dis¬ 
crimination in Section 5. 

2 State discrimination 

In nonlinear quantum mechanics, states that are initially nonorthogonal can evolve to become 
perfectly distinguishable. In this section, we analyze the optimal protocol for distinguishing two 
states of a qubit under the Gross-Pitaevskii and related nonlinearities. 

Specifically, we consider nonlinear dynamics governed by an equation 

= H{tm + K\i^) (2) 

where H{t) is a (time-dependent) Hermitian operator and iF is a nonlinearity of the form 

(x|(a:|V')) = At(Kx|V’)|)(x|V') (3) 

where k: [0,1] —>■ M is a function characterizing the nonlinearity. For example, the Gross-Pitaevskii 
nonlinearity corresponds to k{x) = gx^, where the coefficient g quantifies the strength of the 
nonlinearity. 
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Suppose we are given one of two possible states of a qubit and our goal is to distinguish them as 
quickly as possible. Using the ability to choose the Hermitian driving term H(t), we can orient the 
states however we like, provided we preserve their inner product. Thus, to optimally distinguish the 
given states, we should determine how to orient them on the Bloch sphere so their inner product 
decreases as quickly as possible. 

Consider the pure state with Bloch sphere coordinates {x,y,z), i.e., with density matrix 


1/1 + 2 X — iy\ 

2 \x + iy 1 — z J ' 


(4) 


Since |(0|V’)P = 4^ and |(1|V’)P = 4^, the nonlinear term alone is equivalent to the state- 
dependent Hamiltonian 


0 

In other words, the state evolves according to the equation 


dp 

dt 


1 + 2 x-iy\ /k((4^)^/2) 


k{z) 


X + iy 1 — 2 
0 

ix — y 


IX — y 

0 




where the odd function k{z) : [—1,1] —>■ M is defined by 

R{z) 


(5) 

( 6 ) 

(7) 

( 8 ) 


Thus we find 


d , . 


I^{z) {-y,x, 0 ). 


(9) 


Under these dynamics, states rotate around lines of latitude on the Bloch sphere at a rate depending 
on their latitude. 

Now consider how to optimally orient two states on the Bloch sphere. The rate of change of 
the inner product of Bloch vectors (x+,y+, 2 +) and (x_,y_, 2 _) is 


— (x+x_ +y+y- + 2 + 2 _) = {x+y- -y+x-){K{z+) -k{z-)). (10) 

Suppose the states are separated by a fixed angle a on the Bloch sphere (i.e., angle a/2 in Hilbert 
space) and we aim to rotate them to maximize the rate of decrease of their inner product. Rotations 
about the 2 axis do not affect this rate. Thus, without loss of generality, we can choose the midpoint 
between the two states to lie in the xz plane. We orient the states as shown in Figure 1(a), where 
(j) is the polar angle from the positive 2 axis to the midpoint and 9 is the angle of rotation about 
the midpoint, with 0 = 0 corresponding to the states lying along the line of longitude that passes 
through the x axis. In terms of these parameters, the states have the form 


x± = cos ^ sin (() ± sin ^ cos (/> cos 0 ( 11 ) 

y± = ±sin|sin0 ( 12 ) 

= cos ^ cos 1 ^ + sin ^ sine/)cos 0. (13) 
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Figure 1: (a) Orientation of states on the Bloch sphere with overlap cos f. (b) The flow on the Bloch sphere 
induced by the Gross-Pitaevskii nonlinearity. 


Thus the rate of change of the inner product on the Bloch sphere is 


d 


— cosa = (x+y- - y+x-)(K(z+) - k(z-)) 
= sinQ!sin(/isin0(/t(z_) — k(z+)). 


(14) 

(15) 


Given a specific nonlinearity k, our goal is to choose (p and 9 to minimize (15). Next we perform 
this calculation for several examples. 


2.1 Gross-Pitaevskii nonlinearity 

For the Gross-Pitaevskii nonlinearity, we have k{x) = gx'^, so R{z) = gz and the states evolve as 


d 

dt 


{x,y, 


z) 


gz{-y,x, 0 ). 


(16) 


As depicted in Figure 1(b), we can view this as a flow on the Bloch sphere that pushes states 
along lines of latitude. The rate of rotation varies as a function of latitude, with no rotation at the 
equator or poles and opposite directions of rotation in the northern and southern hemispheres. 

Given a fixed angle between states, we aim to place those states somewhere on the Bloch sphere 
to maximize the rate at which they separate. Equation (15) shows that the rate of change of the 
inner product on the Bloch sphere is 


d 

-- cos a 
dt 


g sin a sin ^ sin^ ip sin 29. 


(17) 


Clearly the rate of decrease of cosa is maximized by choosing p = Trj2 and 9 = Stt/I, giving 


d 

— cos a = —q sm a sm 
dt 


a 

2 • 


(18) 
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Equivalently, the rate of change of the overlap cos ^ (i-e., the magnitude of the inner product in 
Hilbert space) is 


which has the solution 


d „ Id 

— cos 7T = — cos a 

dt 2 4cosfdt 

9 ■ 2 a 

= -Sin 7T 

2 2 


= -|(l-COs2f), 


COS ^ 


cos ^ cosh ^ 


cosh ^ 


- cos 


QQ 


(19) 

( 20 ) 
( 21 ) 


( 22 ) 


where ao is the value of a at t = 0. The states become orthogonal in a time t± such that 
tanh ^ = cos i.e., 


t± 


" 1 
- tanh“ (cos 

9 


2 1 


- ln(cot ^). 
9 


(23) 


In particular, if the initial overlap is cos ^ = 1 — e, the time to distinguish the states is 0(log ^). 
In their optimal orientation, the states have Bloch vectors 


X± = cos J 

i/± = ±^sinf 

.±=±^sinf. 


(24) 

(25) 

(26) 


To keep the states in this orientation, we must apply a rotation about the x axis that keeps the 
y and z components equal. An x rotation (generated by the Hamiltonian gives ^(x,y,z) = 
(z(0, —z,y). Combining this with the effect of the nonlinearity in (16), we have 

^{x,y,z) = {-gyz,gxz - ujz,uy). (27) 

The states remain optimally oriented when 


so we must choose 


o 

II 

1 

(28) 

= gxz - uj{y + z), 

(29) 

(jJ xz . . 

9 y + z 

(30) 

= 2 cos § 

(31) 


(which is given as an explicit function of t by (22)). One can easily verify that, with this choice of 
Lo, equation (27) is satisfied. 

The performance of the optimal procedure using the Gross-Pitaevskii nonlinearity is illustrated 
in Figure 2. As shown in Figure 2(a), if the states are initially close, they separate gradually at 
first, with an accelerating rate of separation, until they become orthogonal. The time for the states 
to become orthogonal (as a function of their initial separation) is plotted in Figure 2(b). 


5 






Figure 2: Optimally distinguishing two states of a qubit using the Gross-Pitaevskii nonlinearity k{x) = gx^. 
(a) Inner product between states as a function of time, as given by equation (22), for initial states with 
oo = 0.1. (b) Time to distinguish states separated by an angle oq on the Bloch sphere, as given by equation 
(23). 


2.2 Logarithmic nonlinearity 

Motivated by a connection to Bose liquids, reference [10] considers the nonlinearity k(x) = 

With this nonlinearity, we have k{z) = s'In so by (15), the inner product on the Bloch sphere 
evolves as 


d ■ ■ j. ■ (l + ^-)(l -^+) 

— cos a = q sin a sin d> sin d In ----. 

dt ^ ^ (l-z_)(l + z+) 


(32) 


Consider the initial states that are optimal for the Gross-Pitaevskii nonlinearity, with (/) = 7r/2 
and 9 = Svr/d. While this choice is suboptimal for the logarithmic nonlinearity, it serves to place an 
upper bound on the time required to distinguish the states. (In fact, numerical calculation shows 
that the optimal states have (p = 7r/2 and 6 ^ Svr/d.) Using these states, we have 


— cos f = g—^ In 


'%/ 2 - 


Sin 


dt 


V2- 


V\/2 + sinf 


sm 2 . 


(33) 


As shown in Figure 3, this is at most a constant times the corresponding quantity under the Gross- 
Pitaevskii nonlinearity, so the performance of the logarithmic nonlinearity is qualitatively similar. 
In particular, for cos ^ = 1 — e we again find that the time to distinguish the states is 0(log ^). 


2.3 General nonlinearities 

So far, we have considered the performance of two specific nonlinearities. However, similar consid¬ 
erations apply for a wide class of nonlinearities of the form (3). 

Since the evolution of states under a nonlinearity k depends only on the function R defined in 
(8), let us evaluate state discrimination in terms of the latter function. Note that we can achieve 
any odd function R for some k. In particular, if we let 


fi{x) X E [0, 

z/(Vl -x2) xE(^,l] 


(34) 
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Figure 3: Rate of change of the overlap for the Gross-Pitaevskii and logarithmic nonlinearities. The solid 
curve is (33) with 5 = 1; the dashed curve is (21) with g = 2. 


for some /r, u: [0, —>■ M, then we find 


k{z) = 



-k{-z) 


z = 0 


ZG [-1,0). 


(35) 


Now suppose R{z) is approximately linear for small z. Specifically, suppose there are constants 
> 0 such that R{z) > gz for all z G [0,5]. For such a nonlinearity, the complexity of distin¬ 
guishing states with overlap 1 — e is 0(^ log ^), just as with the Gross-Pitaevskii and logarithmic 
nonlinearities. 

To see this, a straightforward calculation shows that the states separated by an angle a with 
(j) = 7r/2 and 6 = Stt/A have 


cos f = —sin f) sin f. 
dt 2 ^ VV2 2; 2 


Under the given conditions on k{z), we find 


d a ^ 9 . 

— cos TT <-sm 

dt 2-2 


2 a 
2 


(36) 


(37) 


provided sin ^ < \/25. Comparing with (21), we see that the states separate at least as quickly as 
with the Gross-Pitaevskii nonlinearity provided the overlap is at least Vl — 25^. Similarly to (23), 
a straightforward calculation shows that the time to achieve overlap cos ^ > Vl — 25^ under the 
Gross-Pitaevskii nonlinearity is at most 



cot ^ 
cotf 


(38) 


so the states become distinguishable with constant advantage in time 0{^ log ^). 

In fact, it is not essential for R(z) to be linear about z = 0; similar performance can be achieved 
provided the function grows linearly about any fixed zq G [0,1) (which holds for any function that 
is differentiable and non-constant over some interval). Specifically, suppose there are constants 
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g,A>0 such that |k(zo) ~ + <5)| > ffliJI for all |5| < A. Then, taking cos4> = zq and either 

9 = Svr/d or 0 = 7r/4 (making the choice so that k{zj^) > k{z-)), we find 


Thus by (15), we have 



d ll — Zn 

— cos a < —q\ - - a sin a 

dt - ^ V 2 


(39) 

(40) 


(41) 


provided a is at most some constant (depending on zq and A). Integrating this inequality, we have 
a{t) > e‘^ao where c := g ^/(1 — Zg)/2. Since the angle between the states increases exponentially 
until it reaches at least some constant, states with overlap 1—e become distinguishable with constant 
advantage in time 0{^ log ^). 

Conversely, provided only that R{z) is Lipschitz continuous, this behavior is asymptotically 
optimal. To see this, suppose that R{z) has Lipschitz constant g, i.e., \R{z) — R{z + (5)| < g5. Since 


k+ - z-f < \\{x+,y+,z+) - {x-,y-,x-)f 
= 2(1 — cos a) 

< G? ^ 


(42) 

(43) 

(44) 


we have \i^{z-) — k(2;+)| < ^ga in (15). Therefore 

cosa > — 2g'asinQ!. (45) 

at 

Integrating this inequality gives a{t) < Since the angle increases at most exponentially with 

gt, we require time 0(^ log ^) to distinguish states with overlap 1 — e. 

Note that it is possible to violate this lower bound if the Lipschitz condition does not hold. 
For example, suppose R{z) = sgn( 2 ;) \/\z\ (which is achieved, for example, with k(x) = 0 for 
X G [0, l/-v/2] and k{x) = \/2x‘^ — 1 for x G (l/\/2,1]). Then, taking states with cj) = 7r/2 and 
0 = 3'7r/4, equation (36) gives 


— cos ^ = 
dt 2 


^/ 2 ' 


3/2 


(46) 


While this differential equation does not have a simple closed-form solution, we have the bound 


— cos ^ < —rJ~a sin % 

dt 2 - V 2 


(47) 


for all a G [0,7r], where r := (2^/^y^)“^. Integrating this inequality, we find a{t) > (yTio + rt)'^. 
Therefore any two distinct states become distinguishable in constant time, independent of how 
close they are initially. 









3 Unstructured search 


We now turn our attention to algorithms for unstructured search in nonlinear quantum mechanics. 
In the unstructured search problem, our goal is to search the set [N] := {1, 2,..., A^} for a member 
of the marked set M C [A^], given the ability to determine whether a given element is marked. 
Equivalently, we can consider the problem of deciding whether there is no marked item (i.e., M 
is empty) or a unique marked item (M = {m} for some unknown m G [^])- Using standard 
techniques, an algorithm for this decision problem can be used to find a marked item with only 
logarithmic overhead. 

In the conventional quantum query model, access to the input is provided by a black box that 
determines whether a given element of [A^] is marked (and that can be queried on a superposition 
of elements). Here we work in a continuous-time model [5] where access to the input is provided by 
an oracle Hamiltonian. This Hamiltonian is |m)(m| if M = {m} and is zero if there is no marked 
item. We consider algorithms that use such an oracle Hamiltonian together with an arbitrary 
M-independent driving Hamiltonian H{t) and a given nonlinearity K. The dynamics of such an 
algorithm are governed by the equation 

= {\rn){m\+ (48) 

if M = {m}, or 

+ Kl'ip) (49) 

if there is no marked item. 

We quantify the complexity of such an algorithm by the time required to find the solution. As in 
the standard continuous-time query model, we place no constraints on the norm of H(t). However, 
we consider a fixed-strength nonlinearity K, since an arbitrarily strong nonlinearity could solve 
the search problem arbitrarily fast. Thus our complexities will be expressed as a function of both 
N, the number of items, and g, a parameter characterizing the strength of the nonlinearity. In 
contrast, in the absence of a nonlinearity, the search problem requires II(\/iV) queries even when 
H{t) can be arbitrarily large [5]. 

3.1 Algorithm 

To give an algorithm for unstructured search, we reduce the problem of deciding whether there are 
zero or one marked items to the task of discriminating two possible states of a qubit. We do this 
using the Hadamard test shown in Figure 4. (This approach is similar to that of Abrams and Lloyd 
[3], although we work with the continuous-time query model so we can consider algorithms that 
query the oracle for only a short time.) A straightforward calculation shows that the output of the 
Hadamard test circuit is 

I[|0)(|»> + [/|s)) + |l)(|«)-t/|»»]. (50) 

To produce a single-qubit state, we postselect the second register on the state |s). This postselection 
succeeds with probability 

^[|1 + {s\U\s)\^ + |1 - {s\U\s)\^] = ^ + (51) 
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| 0 ) -{h 


ls) 



Figure 4: The Hadamard test. The first register stores a qubit. The Hadamard gate H acts as H\0) = 
^(|0) + |l))andi/|l) = ^(|0 )-|l))). 


and results in a state 

(l + (s|^|s))|0) + (l-(s|^|s))|l) 
^2(1 + \{s\U\s)\^) 


(52) 


We apply this procedure to the evolution under the oracle Hamiltonian for time ti, namely 
U = if element m is marked, or [/ = / if no item is marked. We choose the uniform 

superposition |s) := initial state. 

If no item is marked, then clearly (s|C/|s) = 1 for any evolution time ti. Thus the postselection 
on |s) always succeeds, and the postselected state of the first qubit is |0). 

On the other hand, if some item m is marked, then we have 


U\s) 


x^m 


x)). 


(53) 


Therefore (s|?7|s) = 1 — ^(1 — e **i), so the success probability (51) is 


l + l(^l^k)P 

2 


= 1 - 


1 

~w~ 


(1 — cosfl) 


1 - 0{tj/N) 


(54) 


(i.e., for ti <C VN, the postselection almost always succeeds). The overlap of the postselected state 
(52) with |0) is 

|l + (^l^|g)l ^ |2iV-l + e--*i| 

^2{l + \{s\U\s)\^) 2^N^-{N-l){l-cosh) 

= ^-^ + 0{tt/N^). (56) 

It remains to distinguish the two possible states. Using the protocol described in Section 2, 
states with overlap 1 — e can be distinguished in time t 2 = 0{^ log ^), where g is the strength of the 
nonlinearity. (We saw in Section 2.3 that such a bound holds not only for the Gross-Pitaevskii and 
logarithmic nonlinearities, but for any nonlinearity of the form (3) where R changes at least linearly 
over some constant-size interval.) From (56), we have e = Q{t\/N‘^). Thus we find an algorithm 
that solves the search problem in total time 


ti+t2 = 0{ti + ^\og^). (57) 

Taking ti = 0(^ \og{Ng)), the total time is 0(^ log(A^g')). 

Of course, if g is very small (in particular, if <C ^^y^) then it may be preferable to eschew 
the nonlinearity and instead use Grover’s algorithm alone. Taking that possibility into account, we 
find an algorithm with complexity 


0(min{^ \og{gN),y/N}) = O 


Vn \ 


.9vW , 1 

logigN) ^ ^ 


(58) 
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This improves the previous upper bound of = 0{min.{^N/g, ^fN}) [9]. For example, 

with g = 0(1), we improve the complexity from 0{y/N) to 0(logA/^). 


3.2 Lower bound 


We now show that the algorithm described above is nearly optimal. We follow the same strategy 
as in the lower bound for the linear case [5]. 

Let \ijj) be the state of the algorithm when there is no marked item and let IV’m) be the state 
when the marked item is m. Consider how the inner product {iplijjm) evolves under the dynamics 
(48) and (49) for an arbitrary driving Hamiltonian H{t). We find 

^(V'l^/’m) = -i('0|m')(m|V’m) + i{K\il)))^\iim) - i('0|(iL|V'm)) (59) 

= -i{\l)\m){m\i)rn) - K(|(x|V’m)|))(V'k)(a^lV'm>. ( 60 ) 

X 

(In particular, observe that the driving Hamiltonian H{t) does not appear, just as in the linear 
case.) Provided |k(x)| < < 7 , we have 


^(k(|(x|V^)|) - K(|(x|V’m)|))(V’k)(x|V’m) < 2c? ^ | (V^| x) (x| | < ‘^9 


by the Cauchy-Schwarz inequality, so 


< |(m|V’)|+ 25 . 


Summing over the N possible marked items, we find 


dt 


■ KV’IV’m)! ^ K"^IV')I) + 25 iv 

m£[N] \me[A''] / 


m 

< Vn + 2gN 


(61) 


(62) 


(63) 

(64) 


where we again used the Cauchy-Schwarz inequality. Integrating for time t using the initial condi¬ 
tion J2m&[N] KV'IV’m)! = N, we find that 

Y, M^m)\>N-tVN{l + 2gVN). (65) 

For the algorithm to succeed with constant probability, the unmarked state must be dis¬ 
tinguishable from each marked state with constant probability, so the final states must satisfy 
\{ijj\'iprn)\ < 1 — 5 for some constant 5 > 0. Therefore 


Y < N{1- 5). ( 66 ) 

m&[N] 


Comparing with (65), we find 


t > 


5Vn 

l + 2g^' 


(67) 


Thus a bounded-error algorithm must have complexity = H(min{2,-y/iV}). This shows 

that the algorithm described in Section 3.1 is optimal up to a logarithmic factor. 
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4 Validity of the Gross-Pitaevskii approximation 


Since the Gross-Pitaevskii model is an approximate description of a fundamentally quantum system, 
upper bounds on the search problem in this model can be interpreted as establishing limitations on 
the effectiveness of the approximation. Using limits on the parallelizability of unstructured search, 
Meyer and Wong argued that this perspective gives a lower bound on how large a Bose-Einstein 
condensate should be for the Gross-Pitaeveskii approximation to hold [9]. They suggested that 
the V-particle Bose-Einstein condensate makes J\f queries in parallel per unit time, so the bound 
= Q{N) for a parallel quantum search algorithm with query depth T [12], together with their 
search upper bound T = 0{^/N/g), implies Af = ^{g). 

While the search algorithm presented in Section 3 has lower complexity, the above argument 
does not apply to that case. In our approach, we query the black box for a short time using only 
linear quantum mechanics and then distinguish the resulting states using the nonlinearity. Since 
the nonlinear portion of the computation makes no access to the oracle, the hardness of parallelizing 
quantum search does not restrict the validity of the Gross-Pitaevskii approximation. 

However, we can use state discrimination to obtain a limitation on the time for which the Gross- 
Pitaevskii approximation remains valid. In this approximation, |'0) is not literally the quantum 
state of the system; rather, it parameterizes the mean-field bosonic state 

where |0) is the vacuum state, aj; is a creation operator for mode x, and Af is the number of particles 
in the system. It is straightforward to understand the distinguishability of these mean-field states in 
terms of their parameterizing vectors: a simple calculation shows that (MF(?/))| M.F((p)) = , 

so the distinguishability of states of the form |MF('0)) is the same as that of In our protocol, 

{ip) is a qubit, so we can represent it using a two-mode condensate and implement the driving term 
H{t) by driving transitions between these two modes. By the Helstrom bound [8], Af = 0(l/e) 
copies are necessary and sufficient to distinguish states with overlap 1 — e with constant advantage 
(in fact, this bound can be achieved even with unentangled measurements [4]). Since the Gross- 
Pitaevskii nonlinearity of strength g can distinguish states with overlap 1 — e in time t = 0{^ log ^), 
we find that the Gross-Pitaevskii approximation can only hold up to time O(^logV). For a 
homogeneous condensate, g = UAf for some constant interaction strength U, so the approximation 
only holds up to time 

While it is appealing to obtain such a limitation from information-theoretic considerations, it 
appears that one can derive a stronger bound by direct analysis of the condensate dynamics [6]. 
This analysis suggests that in fact the approximation is only valid for Ut <C 1/V, so that for a 
homogeneous condensate, it only holds up to time 0{1/Af). 

5 Conclusion 

We have studied the power of nonlinear quantum dynamics to distinguish nonorthogonal states 
and to perform unstructured search. We showed that the Gross-Pitaevskii nonlinearity of strength 
g can be used to distinguish states with overlap 1 — e in time OQ log ^). We applied this result to 
give an algorithm for unstructured search in this model with complexity 0(min{^ \og{gN), y/N}), 
an exponential improvement over previous work. We also showed that many other nonlinearities of 
a related form have similar behavior. Finally, we used our state discrimination protocol to argue 



( 68 ) 
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that for an AA-atom Bose-Einstein condensate, the Gross-Pitaevskii approximation can only hold 
up to time 0(^ logAA). 

We conclude with an open problem regarding nonlinear state discrimination. Our state dis¬ 
crimination analysis was restricted to the case of a single qubit. It is unclear whether states could 
be discriminated more quickly using higher-dimensional systems. Numerical evidence suggests that 
our optimal protocol for distinguishing two states of a qubit with the Gross-Pitaevskii nonlinear¬ 
ity remains optimal if the states are allowed to lie in a higher-dimensional space. Furthermore, 
numerical evidence also suggests that for many other nonlinearities of the form (3), the optimal 
protocol for qubit states cannot be improved by embedding those states in a higher-dimensional 
space. However, we do not have a proof that optimality can be achieved using only two dimensions. 

A notable example where higher dimensions do offer an advantage is provided by the nonlinearity 
with k{x) = x'^ — (as considered in [10]). This nonlinearity has k{z) = 0, so it cannot make 
two states of a qubit more distinguishable. However, by embedding the states in three or more 
dimensions, it becomes possible to decrease their inner product. In this case, numerical evidence 
suggests that there is no advantage to using more than four dimensions. 
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